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PREFACE 


This  paper  presents  the  results  of  a  study  of  Constant 
Absolute  Vorticity  Trajectories  applied  to  isentropic  surfaces. 
The  objectives  of  this  study  were:  first,  to  show  that  Constant 
Absolute  Vorticity  Trajectories  are  theoretically  better  applied 
to  forecasting  future  positions  of  parcels  on  isentropic  surfaces 
than  on  constant  pressure  surfaces  as  they  are  now  applied;  sec¬ 
ondly,  to  show  quantitatively  the  actual  improvement  of  the  tra¬ 
jectories  on  isentropic  surfaces  by  comparing  statistically 
results  of  actual  forecasts  from  both  types  of  charts;  thirdly, 
to  suggest  how  the  improved  technique  can  aid  forecasting,  par¬ 
ticularly  upper  air  prognosis. 

This  work  was  undertaken  as  the  thesis  requirement  for  the 
degree  of  Vaster  of  Science  in  Aerology  at  the  U.  S.  Haval  Post¬ 
graduate  School,  Monterey,  California,  during  the  academic 
year  1952-1953. 

The  author  is  indebted  to  Associate  Professor  George  J. 
Haltiner  of  the  Aerology  Department  for  his  very  valuable  assis¬ 
tance  and  constructive  criticism  during  the  investigation.  He 
also  wishes  to  acknowledge  the  assistance  of  Professor  A.  Boyd 
Mewborn  of  the  Mathematics  Department  and  the  author1 s  wife, 
Pauline,  who  contributed  so  much  help  with  the  laborious  task 
of  decoding,  computing,  and  entering  isentropic  charts. 
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I.  IlTAODUCTIOh 


The  meteorologist  is  often  called  upon  to  forecast  for  periods 
in  excess  of  twenty-four  hours  for  which  simple  extrapolation  is 
generally  insufficient.  One  method  of  forecasting  entails  prepa¬ 
ration  of  prognostic  surface  charts  and  preparation  of  forecasts 
from  these  charts.  As  most  meteorlogists  know,  making  accurate 
prognostic  surface  charts  is  easier  said  than  done.  One  of  the 
important  tools  used  in  estimating  the  position  and  intensity  of 
surface  systems  is  a  prognostic  500  millibar  chart.  At  present, 
much  attention  is  being  paid  to  the  problem  of  upper  air  prognosis, 
and  considerable  information  has  been  published  concerning  fore¬ 
casting  the  500  millibar  surface,  the  most  recent  being  Forecasting 
in  the  Middle  Latitudes  by  Riehl  and  collaborators .  It  is  hoped 
that  this  paper  may  shed  some  light  on  a  facet  of  this  problem  and 
nay  contribute  to  better  prognosis  and  forecasting. 

One  aid  to  prognosis  of  upper  air  charts  is  the  construction 
of  forecast  air  parcel  trajectories  based  on  the  principle  of  con¬ 
servatism  of  the  vertical  component  of  absolute  vorticity,^ +  , 

where  %  is  the  vertical  component  of  relative  vorticity  and  *f 
is  the  vertical  component  of  the  coriolis  force  known  as  the 
coriolis  parameter.  Aossby  to  initiated  the  vorticity  concept 
by  showing  that  under  certain  restrictive  assumptions,  the  follow¬ 
ing  relationship  holds: 


(1) 


The  assumptions  referred  to  above  are  as  follows: 

a.  The  atmosphere  is  barotropic. 

b.  The  atmosphere  is  a  homogeneous,  incompressible  fluid. 

c.  Motion  is  purely  horizontal. 

d.  Friction  forces  are  neglected. 

e.  There  is  no  horizontal  divergence. 

Rossby  £9,  pp  268-289J  then  integrated  his  fundamental  equation 
and  was  able  to  obtain  trajectories  for  particles  which  conserve 
their  vertical  component  of  absolute  vorticity.  This  trajectory 
is  known  as  the  Constant  Vertical  Component  of  Absolute  Vorticity 
Trajectory  and  will  be  referred  to  as  a  CAV  trajectory. 

The  assumptions  necessary  to  derive  the  vorticity  equation 
according  to  Rossby  are  obviously  quite  restrictive.  Starr,  as 
editor  of  the  Journal  of  Meteorology  (June,  1945),  suggested  that 
CAV  trajectories  might  better  be  depicted  on  an  isentropic  chart 
than  on  a  constant  level  chart  as  there  are  no  solenoids  on  an 
isentropic  chart.  Namias  |6,  pp  372-374J  states  that  from  obser¬ 
vations  available  isentropic  surfaces  appear,  to  the  first  approxi¬ 
mation,  to  be  substantial  surfaces,  i.e.,  surfaces  that  contain  the 
same  air  particles  from  day  to  day.  Influences  of  non-adiabatic 
nature  are  appreciable  over  long  periods  of  time,  but  are  usually 
insufficient  to  disrupt  the  fundamental  isentropic  flow  from  one 
day  to  the  next.  The  derivation  of  the  vorticity  equation  for  flow 
on  an  isentropic  surface  requires  fewer  assumptions.  Such  restrictions 
as  purely  horizontal  flow,  a  barotropic  atmosphere,  and  a  homogeneous 
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incompressible  atmosphere  may  be  removed.  Now,  an  isentropic 
surface  is  a  surface  of  constant  potential  temperature.  In 
meteorlogical  processes,  potential  temperature  is  one  of  the  more 
conservative  elements  by  which  parcels  of  air  nay  be  identified 
from  time  to  time.  Since  the  atmosphere  is  normally  stable, 
potential  temperature  increases  with  height  and  the  atmosphere 
may  be  considered  to  consist  of  an  infinite  number  of  isentropic 
surfaces.  Sir  Napier  Shaw  £l0,  p  263]  first  suggested  that  charts 
of  isentropic  surfaces  should  be  drawn  as  the  motion  of  air  is  best 
resolved  on  isentropic  surfaces.  It  will  be  shown  in  this  paper 
that  the  prognostic  trajectory  based  on  the  conservatism  of  the 
vertical  component  of  absolute  vorticity  is  theoretically  and 
practically  better  suited  to  use  on  isentropic  charts  than  on 
constant  pressure  charts  as  it  is  now  applied. 
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II.  THEORETICAL  INVESTIGATION 


The  motion  of  particles  on  a  substantial  surface  can  be 
described  with  Lagrangian  methods.  The  form  of  the  vorticity 
equation  for  flow  on  an  isentropic  surface  derived  here  using  the 
Lagrangian  method  of  attack  is  due  to  Ealtiner*'.  LetX,ij#Z  and-fc 
be  normal  Cartesian  coordinates  and  let^X^X,  Z?  and  y  be 
Lagrangian  coordinates  where  ~X  and  jC  are  coordinates  of  the 
horizontal  projection  of  the  isentropic  surface,  ~  4  and  ^ 
refers  to  time  for  the  isentropic  surface.  Then 
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It  follows  from  the  above  equations  that 
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the  stream  function  for  an  isentropic  surface. 


In  a  similar  way  we  find 
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Further,  since  the  motion  is  assumed  to  be  isentropic, 
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Substituting 

in 

the  equations  of  notion 
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the  expressions 
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,  we  obtain 
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Differentiating  (7)  with  respect  to  T  and  (8)  with  respect  to  X 
and  subtracting,  we  obtain 
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and  introducing  vector  notation  we  obtain 


-  (£[£]  %)  ~fa  V  v= +iviev  +az 

where  the  operator  ^  =  £j  ^ _ ,  n  o> 

^  c>X  ^ 

Finally,  we  may  write  this  as 

ikfy+t)+(r0+f)*z»v=o  (S) 

since"}"  — t  for  any  particle.  Now,  assume  the  divergence  tern  -V=0. 

£& 

Then  we  may  write  the  vorticity  equation  in  the  familiar  form 


While  this  form  is  sirailar  to  the  form  derived  by  Ilossby,  there 
are  some  fundamental  differences.  Parcel  motion  is  not  limited 
to  the  horizontal,  but  is  three  dimensional.  Moreover,  a  barotropic 
and  incompressible  atmosphere  is  not  required. 
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III.  TECHNIQUE  OF  INVESTIGATION 


The  relative  vorticity  may  be  written  in  a  form  similar  to 
the  usual  expression  for  the  vertical  component  of  vorticity  in 
natural  coordinates 

y  _  v  _  av 

OH 

where  P  and  fsj  are  measured  on  the  isentropic  map.  If  we  select 
a  parcel  in  a  broad  uniform  flow  or  in  the  axis  of  a  jet  stream, 

dV 

Further,  if  the  parcel  is  at  an  inflection  point  in  the  flow,  hL  —  O. 

*5 

Under  these  conditions  the  CAV  trajectory  can  be  computed.  It 
will  be  noted  that  defined  above  is  not  exactly  equal  to  the 
projected  radius  of  curvature  of  the  streamline  on  the  isentropic 
surface.  However,  the  approximation  involved  here  is  normally 
less  than  those  made  in  the  actual  computation  of  CAV  trajectories 
as  described  by  Fultz  [4,  P  13]  and  revised  by  Wobus  [12]  .  \vhile 
approximations  have  to  be  made,  it  must  be  remembered  that  the  true 
worth  of  a  technique  lies  in  the  statistically  proven  results 
obtained. 

The  simplest  method  of  obtaining  a  CAV  trajectory  is  through 
the  use  of  a  mechanical  device  such  as  the  so-called  ’’Wiggle-Wagon” 
devised  by  Wobus  £l2j  and  currently  in  use  in  the  WBAN  Analysis 
Center,  Washington,  D.  C.  Original  studies  of  the  CAV  trajectory 
assumed  a  plane  earth.  This  assumption  entailed  errors  due  to 
distortion  when  compared  to  trajectories  along  a  spherical  earth. 

The  Wobus  ’’Wiggle-Wagon”  corrects  this  defect.  To  aid  meteorolo¬ 
gists  who  do  not  have  access  to  the  ”Wiggle-Wagon”,  the  U.  S.  Navy 
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Bureau  of  Aeronautics  Project  AROWA  has  published  a  table 

for  computing  CAV  trajectories  made  fron  trajectories  traced  by 
the  ’Wiggle-Wagon” .  This  table  will  be  used  to  compute  CAV  tra¬ 
jectories  in  this  paper.  For  application,  it  is  necessary  to  know 
speed,  direction,  and  the  latitude  of  the  parcel,  all  at  an  inflection 
point  in  the  flow. 

The  next  step  in  the  investigation  was  to  compute  trajectories 
on  isentropic  and  500  millibar  charts  and  test  the  results  statisti¬ 
cally.  A  rather  complete  compilation  of  upper  air  and  surface  data 
are  available  in  the  Northern  Hemisphere  Surface  and  500  Itillibar 
Charts  Series  published  by  the  U.  S.  Weather  Bureau.  These  remarkable 
publications  not  only  have  a  daily  series  of  logically  analysed  surface 
and  500  millibar  charts,  but  also  contain  a  compilation  of  all  data 
used  to  plot  these  charts.  Data  were  taken  from  the  period  of 
January-February  1949.  From  this  data,  a  set  of  thirty  isentropic 
charts  was  plotted.  The  303K  isentropic  surface  was  selected  because: 
(1)  this  surface  was  often  close  to  the  so-called  "level  of  non¬ 
divergence”  described  by  Bjerknes  ;  (2)  this  surface  often  con¬ 
tained  much  of  the  500  millibar  wind  speed  maximum  axis  (jet); 

(3)  the  303K  isentropic  surface  was  usually  sufficiently  far  from 
the  surface  to  remove  effects  of  surface  friction  and  turbulence. 

It  was  desired  to  select  an  isentropic  surface  containing  much  of 
the  500  millibar  isotach  maximum  axis  so  computation  of  CAV  tra¬ 
jectories  could  aid  in  the  prognosis  of  the  500  millibar  surface. 
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Inasnuch  as  the  derivation  of  CAV  trajectories  on  an  isentropic 
surface  assumes  the  divergence  tern  to  be  aero,  the  selection  of 
an  isentropic  surface  close  to  the  level  of  no  horizontal  divergence 
nay  be  desirable.  The  effect  of  the  divergence  tern  is  discussed  in 
Chapter  IV. 

After  the  heights  of  the  30311  isentropic  surface  were  plotted, 
all  available  wind  data  pertinent  to  the  surface  were  plotted.  Where 
wind  observations  were  available  from  stations  other  than  radiosonde 
stations,  an  interpolated  value  of  the  height  of  the  isentropic  surface 
and  the  wind  nearest  that  height  were  entered.  Because  of  the  fairly 
dense  network  of  upper  air  word  observing  stations  in  the  United 
States,  sufficiently  accurate  streamline  analysis  was  possible  with¬ 
out  the  time-consuming  computations  of  stream  function  values. 

Streamline  analyses  were  made  for  0300Z  each  day  in  the  period 
from  6  January,  1S49,  to  22  January,  1949,  and  from  10  February,  1949, 
to  20  February,  1949.  For  each  chart,  a  number  of  CAV  trajectories 
were  computed.  The  following  conditions  had  to  be  met  before  choosing 
the  initial  point:  (1)  the  point  had  to  be  in  or  very  near  an  inflection 
point  in  the  flow  (undergoing  no  effects  of  curvature);  (2)  the  point 
had  to  be  in  or  near  the  axis  of  the  isotach  maximum  (jet  stream); 

(3)  the  point  must  be  at  an  actual  point  of  wind  observation  to  obtain 
correct  initial  wind  speed  and  direction.  This  procedure  netted  from 
one  to  four  points  suitable  for  computation  of  CAV  trajectories  per 
chart. 
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For  comparison  with  each  CAV  trajectory  plotted  on  the 
isen tropic  chart,  a  companion  CAV  trajectory  was  computed  on  the 
500  millibar  chart  of  the  sane  time.  It  was  not  too  difficult 
to  find  a  companion  inflection  point  on  the  500  millibar  chart  as 
the  isentropic  surface  was  usually  close  to  the  500  millibar  surface. 
After  the  CAV  trajectories  were  constructed  on  each  chart,  forecast 
wind  velocities  were  determined  by  moving  the  parcel  along  the  tra¬ 
jectory  a  distance  equal  to  the  initial  wind  speed  multiplied  by 
the  number  of  hours  to  forecast  time.  The  forecast  times  selected 
were  twenty-four,  forty-eight,  and  seventy-two  hours.  These  fore¬ 
cast  winds  were  then  plotted  on  the  verifying  charts  and  verified. 

On  the  isentropic  chart,  verification  depended  on  the  forecast 
position  falling  within  the  observational  network.  If  the  forecast 
position  fell  within  the  network,  but  not  on  a  particular  station, 
the  forecast  winds  were  verified  by  linear  interpolation  between 
nearby  stations.  However,  care  was  taken  to  use  only  stations 
where  linear  interpolation  would  be  reasonable.  On  the  500  millibar 
chart,  forecast  winds  were  verified  in  a  similar  manner  if  the  fore¬ 
cast  position  fell  within  the  observational  network  or  on  a  reporting 
station.  If  the  forecast  position  on  the  500  millibar  chart  fell  on 
neither  of  these  places,  then  contours  were  assumed  to  be  streamlines, 
and  gradient  winds  were  measured  to  verify.  Since  verification  on 
isentropic  charts  depended  on  forecast  positions  falling  in  the  North 
American  network  or  on  a  particular  station,  the  number  of  possible 
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verifications  fell  off  drastically  after  forty-eight  hours.  However, 
with  hemispherical  500  millibar  charts  available,  forecasts  made  on 
the  500  millibar  chart  could  always  be  verified.  All  500  millibar 
CAV  trajectories  were  carried  out  to  seventy-two  hours. 

In  verifying  forecast  winds,  observed  wind  direction  error 
was  arbitrarily  selected  negative  if  the  observed  wind  direction  was 
to  the  right  (clockwise)  of  the  forecast  wind  direction  and  positive 
if  the  observed  wind  direction  was  to  the  left  (counter-clockwise) 
of  the  forecast  wind  direction.  Similarily,  wind  speed  error  was 
selected  negative  if  observed  wind  speed  was  less  than  forecast  and 
positive  if  the  wind  speed  was  greater  than  forecast.  All  forecast 
winds  and  verifications  were  tabulated  in  a  fora  suitable  for 
statistical  testing. 

Histograms  of  wind  speed  error  and  wind  direction  error  were 
prepared  for  twenty-four  and  forty-eight  hour  forecast  times  for  both 
the  isentropic  and  500  millibar  charts.  These  histograms  are  depicted 
as  Figures  1  through  4.  Wind  speed  errors  were  grouped  into  cells 
each  having  a  cell  interval  of  10  knots  error.  The  cell  midpoints 
are  fori  0,^10, ±.ZO  etc.,  knots  error.  Since  the  cell  boundaries 
are  exactly  halfway  between  the  cell  midpoints,  some  of  the  obser¬ 
vations  fell  on  the  cell  boundaries.  In  this  case,  it  was  deemed 
best  to  divide  such  observations  between  each  adjacent  cell.  Less 
trouble  was  encountered  in  grouping  wind  direction  error.  /Ill  wind 
direction  errors  were  in  increments  of  10  degrees  and  values  of O,±!0 
±ZO,  etc .  degrees  of  error  were  taken  as  cell  midpoints • 
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After  the  data  were  grouped  in  histograms,  normal  curves  were 
fitted  in  accordance  with  the  method  described  by  Hoel  £5,  pp  191-194} 
After  obtaining  these  fitted  normal  curves,^  tests  were  performed  to 


determine  if  each  sample  of  data  could  reasonably  be  from  a  normal 
population.  In  all  cases  the  sample  distributions  could  reasonably 
be  said  to  be  from  normal  populations,  although  in  two  cases  the 


value . 


test  were  raet  by  grouping  cells  with  too  few  frequencies  together. 

After  showing  normality  of  the  data,  the  distributions  of  data  from 
the  isentropic  charts  were  compared  to  similar  distributions  from 
the  500  millibar  charts  by  the  use  of  the  UFM  test.  The  f,F"  test 
was  made  to  determine  if  compared  data  distributions  from  the 
isentropic  and  500  millibar  charts  could  statistically  be  from  the 
same  parent  normal  population.  The  purpose  here  is  to  show  that  a 
given  distribution  of  data  from  the  isentropic  charts  is  not  from  the 
same  normal  population  as  the  similar  data  distribution  from  the  500 
millibar  charts.  If  there  is  a  significant  difference,  then  inspection 
and  intuition  will  reveal  which  of  the  two  types  of  charts  are  better 
suited  to  forecasting  by  means  of  CAV  trajectories. 
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IV.  RESULTS 


It  is  well  to  look  first  at  the  data  in  "everyday”  terns  to 
find  what  kind  of  results  a  forecaster  nay  expect.  If  it  is  arbi¬ 
trarily  decided  that  a  wind  direction  forecast  will  verify  as  a  "hit" 
if  it  is  within  20  degrees  of  observed  direction,  then  we  find  that 
96.4%  of  the  twenty-four  hour  wind  direction  forecasts  from  GAV  tra¬ 
jectories  on  isentropic  charts  verified  as  "hits".  By  comparison, 

76.3 %  of  similar  forecasts  frora  500  millibar  charts  verified  as  "hits". 
With  forty-eight  hour  wind  direction  forecasts,  81,0  frora  isentropic 
charts  verified  as  "hits"  and  41*  frora  500  millibar  charts. 

If  it  is  arbitrarily  decided  that  15  knots  of  wind  speed  error 
is  not  too  much  to  verify  as  a  "hit",  we  find  that  76%  of  twenty-four 
hour  wind  speed  forecasts  from  trajectories  on  isentropic  charts 
verified  as  "hits".  Of  the  data  frora  the  500  millibar  charts,  42.4,o 
of  twenty-four  hour  wind  speed  forecasts  verified  as  "hits".  Similarly, 
forty-eight  hour  wind  speed  forecasts  verified  78.5/o  from  isentropic 
chart  data  and  39,*  from  500  millibar  chart  data.  So,  percentagewise, 
forecasting  wind  vectors  by  CAV  trajectories  on  isentropic  charts 
was  superior  to  forecasting  wind  vectors  by  GAV  trajectories  on  500 
millibar  charts  during  the  period  of  this  study. 

Figures  1  through  4  are  histograms  of  error  between  observed 
and  forecast  wind  speed  and  direction.  Applying  standard  statistical 
methods,  normal  curves  were  fitted  to  the  data.  The  test  was 
applied  to  determine  if  the  distributions  of  error  were  likely  to  be 
from  a  normal  population.  Only  the  twenty-four  and  forty-eight  horn" 
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errors  were  compared  and  tested  as  the  seventy- two  hour  isentropic 
sample  was  too  small.  In  each  case  tested,  the  distribution  was 
shown  to  be  reasonably  from  a  normal  population.  In  cases  of  data 
from  isentropic  charts,  there  was  little  doubt  that  the  distributions 
were  from  normal  populations.  But  in  cases  from  the  500  millibar 
charts,  the  values  of^(  were  close  to  the  critical  values  and, 
therefore,  there  is  some  doubt  that  they  may  be  from  a  normal  popu¬ 
lation.  However,  since  the  values  of  did  fall  below  the  critical 
values,  we  will  treat  such  data  as  being  from  normal  populations. 

The  "F"  test [5,  pp  152-154]  is  a  statistical  test  for  comparing 
variances  (s')  of  two  samples  of  data  to  determine  if  the  two  samples 
can  reasonably  be  from  the  same  normal  population.  The  basic  require¬ 
ment  of  the  nFM  test  is  that  the  data  samples  must  each  be  from  normal 
populations ,  although  not  necessarily  the  same  normal  population.  This 
is  why  it  was  necessary  to  show  normality  of  the  data  samples  above. 

To  apply  the  "Fn  test,  we  first  postulate  that  there  is  no  significant 
difference  between  samples  tested,  i.e.,  S,X  =■  5 J”  .  This  is  the  null 
hypothesis.  If  the  value  of  F  computed  falls  below  a  critical  value, 
the  null  hypothesis  is  not  denied  and  there  is  no  statistical  difference 
between  the  two  samples.  If  the  value  of  F  computed  is  higher  than  the 
critical  value,  the  null  hypothesis  is  denied,  and  there  is  a  statisti¬ 
cally  significant  difference  between  the  two  samples  and  they  are  not 
reasonably  from  the  same  parent  normal  population.  It  is  this  last 
condition  we  wish  to  obtain  here.  The  results  of  applying  the  ,TFn 
test  are  as  follows: 
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1.  Twenty- f our  Hour  Wind  Direction  Error.  The  variance  of 
data  from  the  isentropic  charts  5,  L-=  M4  .  The  variance 
of  data  frori  the  500  millibar  charts  5^=:  6.37  .  Then 

F=  ST/S,l= 

at  a  95/c  level  of  belief  is  4.90.  Therefore,  the 
samples  tested  are  significantly  different  and  the  null 
hypothesis  is  denied. 

2.  Forty-eight  Hour  >/ind  Direction  Error.  The  variance  of 

date  from  the  isentropic  charts  5^~~  3.  1 5  .  The 

variance  of  data  from  the  500  millibar  charts  5^—30.12.  . 
Then 

F~  —  5.7  ST=  r  IS. 

at  a  95,0  level  of  belief  is  3.75.  This  is  highly 
significant  and  the  null  hypothesis  is  certainly  denied. 

3.  Twenty-four  Hour  hind  Speed  Error.  The  variance  of  data 
from  the  isentropic  charts  is  155.7.  The  variance  of  data 
from  the  500  millibar  charts  is  413.1.  Then 

F  -  2XS. 

F<_  at  a  95,0  level  of  belief  is  5.19  and,  therefore,  the 
null  hypothesis  is  not  denied  and  there  is  no  significant 
difference  between  the  samples  cf  data. 

4.  Forty-eight  Hour  Vind  Speed  Error.  The  variance  of  data 
frori  the  isentropic  charts  is  127.3.  The  variance  of  data 

from  the  500  millibar  charts  is  515.5.  Then 

F  =  4.05. 

At  a  90, i  level  of  belief  the  value  of  4.00  ,  and  at 
this  level  the  difference  between  the  samples  becomes 
significant. 


(15) 


Significant  results  were  obtained  in  the  cases  of  wind  direction 
error  and  in  one  case  of  wind  speed  error.  In  each  of  the  cases 
where  significance  was  obtained,  an  intuitive  examination  of  the 
data  shows  that  it  is  the  data  from  trajectories  constructed  on 
isentropic  charts  that  are  better. 

It  is  gratifying  to  note  that  the  greatest  improvement  was 
shown  for  i/ind  direction  forecasting  for  a  forty-eight  hour  period. 

In  nuch  operational  forecasting,  wind  direction  is  usually  deemed 
more  important  than  wind  speed  for  prognostic  work.  The  proper 
placing  of  troughs  and  ridges  and  associated  weather  patterns 
depends  greatly  on  correct  wind  direction  forecasts.  Therefore, 
from  the  results  obtained,  the  CAV  trajectory  method  will  show 
greatest  improvement  over  the  present  system  of  taking  CAV  tra¬ 
jectories  on  constant  pressure  charts  if  used  on  isentropic  charts 
to  forecast  wind  direction,  at  least  up  to  forty-eight  hours. 

A  recent  study  of  CAV  trajectories  along  a  psuedo-600  millibar 
chart  was  published  by  Bruch  DO  .  Bruch’s  results  on  wind  direction 
error  are  quite  similar  to  results  obtained  herein  for  the  500  millibar 
chart.  Bruch  did  not  verify  wind  speeds. 

It  is  realized  that  the  isentropic  chart  leaves  much  to  be  desired 
as  an  operational  weather  chart.  At  present  the  chart  is  of  course 
laborious  to  plot  since  no  isentropic  data  are  transmitted  on  tele¬ 
type  circuits;  however,  if  such  data  were  ti'ansnitted  on  teletype 
circuits,  streamline  analysis  would  take  little  longer  to  prepare 
than  the  present  500  millibar  chart,  /mother  difficulty  is  that  the 
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Since 


nover. ent  of  the  isentropic  surface  is  rot  easy  to  predict, 
the  CAV  trajectory  developed  herein  is  a  tv.ro  dimensional  projection 
of  a  three  dimensional  path,  the  height  of  the  parcel  in  the  future 
is  in  doubt.  Of  course,  if  the  isentropic  sxirface  did  not  r.ovc,  the 
future  position  of  a  parcel  would  be  certain.  But,  this  is  not 
generally  the  case.  A  future  study  to  correlate  this  study  with  the 
prognosis  of  isentropic  surfaces  could  be  to  advantage,  honetheless, 
the  computation  of  CAV  trajectories  on  an  isentropic  surface  can  be 
of  value  to  prognosis  of  a  nearby  constant  pressure  surface.  The 
positions  of  prognostic  troughs  and  ridges  of  the  isentropic  flow 
can  be  placed  on  a  nearby  constant  pressxu-e  surface. 

From  study  of  the  CAV  trajectories  on  the  isentropic  charts  of 
this  paper,  it  qualitatively  appears  that  parcels  imbedded  in  north¬ 
westerly  flow  descend  as  they  move  southeastward,  whereas  parcels 
imbedded  in  southwesterly  flow  rise  as  they  move  northeastward. 

In  the  final  form  of  the  vorticity  equation  there  is  a  diver¬ 
gence  tern  that  is  assumed  to  be  zero.  This  assumption  is  not  always 
fulfilled  in  the  atmosphere.  The  effect  of  this  divergence  term  may 
be  important.  For  illustration,  let  the  divergence  tern 
be  constant  in  equation  (9)  and  integrating  with  respect  to  -fc  ,  we 
obtain 
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Let  *£rrO  at  the  initial  inflection  point,  where  alsoj^=o  . 


(Tc,+*)  =  c  = 


and 
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As  can  be  seen,  positive  values  of  the  divergence  tern  in  a  north¬ 
westerly  flow  will  cause  a  negative  contribution  to  relative 
vorticity  of  a  parcel  causing  it  to  go  farther  south  than  under  a 
CAV  trajectory.  Positive  values  of  the  term  in  southwesterly  flow 
will  have  a  similar  effect  on  the  relative  vorticity  of  a  parcel 
causing  it  to  not  reach  the  maximum  latitude  indicated  by  a  CAV 
trajectory.  Negative  values  of  the  divergence  tern  will  cause 
positive  contributions  to  the  relative  vorticity  of  a  parcel.  In 
northwesterly  flow  this  will  cause  a  parcel  to  curve  more  cycloni- 
cally  than  would  be  indicated  by  a  CAV  trajectory  and  to  curve  less 
cyclonically  in  southwesterly  flow.  Fultz  jj4,  pp  32-109J  mentions 
that  in  a  lower  level  (10,000  foot  level)  the  average  type  of 
deviation  from  CAV  trajectories  is  such  that  could  be  caused  by 
horizontal  velocity  divergence  in  northwesterly  flow  and  convergence 
in  southwesterly  flow.  Mow, 


Oi&X  c>Y/  . 


Thus,  the  divergence  term  is  composed  of  two  parts,  horizontal 


velocity  divergence  and  a  shear  tern 


.  Fleagle  and  Panofsky  £*7^ 
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have  given  several  values  for  horizontal  velocity  divergence, 
with  a  range  of  about ^  5  x  lCT^  per  second.  The  values  of  the 
shear  tern  could  be  computed  at  initial  inflection  points  for 
thirty-four  of  the  CAV  trajectories  on  isentropic  charts  of  this 
study.  In  fourteen  cases  the  terra  was  less  than  10“  ^  per  second. 

In  eleven  cases  the  value  of  the  shear  term  was  approximately  10“®, 
and  in  nine  cases,  approximately  -10“®  per  second.  Usually,  the 
shear  term  was  positive  or  negative  if  the  flow  was  up  or  down  the 
isentropic  surface  respectively.  From  this  we  see  that  in  about 
60fa  of  the  cases  studied,  the  shear  tern  was  of  the  sane  order  of 
magnitude  as  the  horizontal  velocity  divergence.  In  these  cases 
the  shear  tern  either  canceled  or  doubled  the  horizontal  velocity 
divergence  term .  Now,  if  the  assumption  ^26>'  ^  °  were  the 

most  important  assumption  in  the  development  of  CAV  trajectories,  we 
roight  expect  forecasting  of  future  positions  of  parcels  by  CAV  tra¬ 
jectories  on  isentropic  surfaces  to  give,  in  some  cases,  greater 
error  than  from  applying  CAV  trajectories  to  constant  pressure  sur¬ 
faces.  This  was  not  observed  in  the  data  collected  in  this  study. 
This  nay  sttggest  that  the  assumption  of  no  horizontal  divergence  is 
not  too  restrictive,  and  that  possibly  one  of  the  other  assumptions 
made  by  Rossby  is  more  restrictive,  e.g.,  assumption  of  purely 
horizontal  notion. 

If  we  take  a  large  constant  value  for  horizontal  velocity 
divergence,  say  6  x  10“^  per  second,  take  a  like  value  for  the 
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shear  terra,  assume  a  time  period  of  twelve  hours,  and  use  these 
values  in  equation  (11),  we  would  find  that  the  value  of  at, 
say,  45°  is  reduced  by  33, o  or  to  the  value  of-f  at  latitude  33°. 
But,  this  is  an  extreme  case,  and  the  effects  of  the  divergence 
term  are  usually  much  less.  However,  since  the  assumption  of 
V,ev=o  does  lead  to  error,  it  seems  that  some  future  re¬ 
search  could  be  profitably  spent  introducing  quantitative  values 
for  the  divergence  term  and  computation  of  a  trajectory  incorpor¬ 
ating  this  correction  to  the  GAV  trajectory  for  use  in  more 
accurate  forecasting. 
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FREQUENCY 


ERROR 


Distribution  of  twenty-four  hour  wind  direction 
error  from  (a)  isentropic  charts  and  (b)  from 
500  millibar  charts 
Zrror:  tens  of  degrees 

Figure  1. 
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FREQUENCY 


Distribution  of  forty-eight  hour  wind  direction 
error  from  (a)  is en tropic  charts  and  (b)  from 
500  millibar  charts 
Error:  tens  of  degrees 

Figure  2. 
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FREQUENCY 


V 


Distribution  of  twenty-four  hour  wind  speed 
error  from  (a)  isentropic  charts  and  (b) 
from  500  millibar  charts 
Error :  knots 

Figure  3. 
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FREQUENCY 


20 


Distribution  of  forty-eight  hour  wind  speed 
error  from  (a)  isentropic  charts  and  (b) 
from  500  millibar  charts 
Error:  knots 


Figure  4. 
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